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In 1892, J. A. Allen named Perognathus merriami on the 
basis of 17 silky pocket mice from the vicinity of Brownsville, 
Texas. He compared this material with specimens of P. flavus 
from El Paso and concluded that they represented distinct 
species. He also figured six skulls of P. merriami and three 
of P. flavus. 

Allen used the following characters to distinguish P. 
merriami from P. flavus: “. . . brighter and more yellowish in 
coloration, the sides being strongly yellowish or golden in- 
stead of pale cinnamon. . . . generally brighter yellowish 
color, particularly along the sides. . . . the skull is much 
larger and disproportionately broader in proportion to the 
length. The mastoids are shorter and less developed, leaving 
a much broader intermastoid area, with the interparietal much 
broader than long, instead of nearly square as in P. flavus." 

W. H. Osgood (1900), in his revision of the genus Perog- 
nathus, named Perognathus merriami gilvus based on three 
specimens from west Texas and four from Eddy (= Carls- 
bad), New Mexico. Osgood recognized the difficulty of sepa- 
rating his specimens from P. flavus: "This subspecies com- 
bines to some extent the characters of flavus and merriami” 

I have studied more than 250 specimens of P. merriami 
and P. flavus, including the types of P. merriami, P. merriami 
gilvus, P. flavus bimaculatus, P. f. fuliginosus, P. f. hopiensis, 
P. f. mexicanus, P. f. piperi, and P. f. sonoriensis. I report here 
the results of a discriminant function analysis which leads me 
to believe that P. merriami and P. flavus are conspecific. 
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Fic. 1. Map showing collecting localities of Perognathus merriami 
Wl, Perognathus flavus @, and Perognathus merriami gilvus A. ® indi- 
cates locality of supposed sympatry of P. flavus and P. m. gilvus. 


METHODS 


Figure 1 shows the specimen localities listed in Table 1. 
All specimens are study skins, skulls, or both. A total of 24 
characters was used in the discriminant function analysis; 
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TABLE 1. 


Sample Locality 


Samples of Perognathus examined 


Used in 
Multivariate 
Analysis 


Perognathus merriami gilvus 


Mexico, Coahuila, 15SE Langtry, Texas 
New Mexico, Carlsbad 


oe 


New Mexico, Roswell, 40W 


New Mexico, Tucumcari, 25W 


Texas, Presidio Co. 
Texas, Langtry 
Texas, Comstock 
Texas, Stanton 

Texas, Washburn 
Texas, Big Springs 
Texas, Fort Lancaster 
Texas, Marathon, 15S 
Texas, Marathon 
Texas, Boquillas 
Texas, Juno, 20E 
Texas, Rock Spring 
Texas, Japonica, 15W 
Texas, Howard, 5S 
Texas, Monahans 
Texas, Comstock, 10S 
Texas, Van Horn 
Texas, Samuels 
Texas, Eagle Pass 
Texas, Mobeetie 
Texas, Kermit, E 
Texas, Magic City, IN 


Ë = O = Pee 


Subtotals 9 


Perognathus merriami merriami 
Mexico, Coahuila, Monclova 2 
Mexico, Coahuila, Sabinas 
Mexico, Coahuila, Saltillo 
Mexico, Nuevo Leon, Linares 1 
Mexico, Nuevo Leon, Doctor Cos 1 
Mexico, Nuevo Leon, Aldama 


Mexico, Tamanlipas, C. Victoria 
Mexico, Tamaulipas, Hidalgo 
Mexico, Tamaulipas, Mier 

Mexico, Tamaulipas, Matamoros, 10S 


Nw QL > 


Mexico, Tamaulipas, Reynosa 
Mexico, Tamaulipas, Alta Mira 
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Additional 


LO = = oO 


bo r2 op 
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TABLE I, 


Sample Locality 


(Continued ) 


Mexico, Tamaulipas, Nuevo Laredo 


Mexico, Tamaulipas, Camargo 


Mexico, Tamaulipas, Soto La Marina 


Texas, Brownsville 
Texas, Blocker Rancho 
Texas, Rio Grande City 
Texas, Padre Island 
Texas, Santa Rosa 
Texas, Kerrville 

Texas, Brownsville, 20E 
Texas, Mason 


Perognathus flavus flavus 
Mexico, Chihuahua, Gallego 
Mexico, Chihuahua, Chihuahua 
Mexico, Chihuahua, Escalon 
New Mexico, Deming 
New Mexico, Chico Springs 
New Mexico, Dog Spring 
New Mexico, Carrizozo 
New Mexico, Tularosa 
New Mexico, Santa Rosa 
New Mexico, Ribera 
New Mexico, Mesa Jumanes 
New Mexico, Manzano Mtns. 
New Mexico, Grants 
New Mexico, Carlsbad 
New Mexico, Rio Puerco 
New Mexico, Quemado, 10SW 
New Mexico, Glenwood 
New Mexico, Laguna 
New Mexico, Hachita 
New Mexico, Playas Valley 


Subtotals 


Used in 
Multivariate 
Analysis 

2 
1 
2 
29 


= t° 


New Mexico, Gila National Forest, Diamond Creek 1 


New Mexico, Fairview 
New Mexico, Reserve, 5NE 
New Mexico, Vaughn, 5W 
Texas, El Paso 

Texas, Sierra Blanca 
Texas, Valentine 

Texas, Alpine 


Subtotals 


Additional 
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TABLE 1. (Continued) 


Used in 


Multivariate 
Sample Location Analysis 

Perognathus flavus medius 

Mexico, San Luis Potosi, Hacienda La Parada 1 

Mexico, San Luis Potosi, Jesus Maria 1 
Perognathus flavus bunkeri 

Colorado, Loveland 1 
Perognathus flavus hopiensis (type) 

Arizona, Oraibi 1 


Perognathus flavus mexicanus (type) 
Mexico, Distrito Federal, Tlalpam 
Perognathus flavus piperi (type) 
Wyoming, Newcastle 1 
Perognathus flavus sonoriensis (type) 
Mexico, Sonora, Costa Rica Ranch 
Perognathus flavus fuliginosus (type) 
Arizona, San Francisco Mtn. 
Perognathus flavus bimaculatis (type) 
Arizona, Fort Whipple 


Subtotals 2 
"Totals 136 


179 


Additional 


L2 r3 


nine were skull characters measured to the nearest .05 mm 
with dial calipers, eight were skin measurements or scores, and 
seven were ratios (Table 2). The characters used were as 


follows: 


Occipitonasal length.—Medial distance from the tip of the 
nasals to the posteriormost margins of the condyles. 
Frontonasal length.—Medial distance from the tip of the 


nasals to the frontoparietal suture. 


Nasal length.—Medial distance from the tip of the nasals 


to the nasal-frontal suture. 


Least interorbital constriction.—The least distance across 


the skull between the orbits. 


Mastoid breadth.—Greatest width of the skull across the 


mastoid bones. 


Bullar length.—Greatest length of the auditory bulla. 
Interparietal length.—Greatest length of interparietal bone. 
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Fic. 2, Five standards used for scoring skins for amount of dorsal 
black and degree of development of postauricular patches. 


Interparietal width.—Greatest width of interparietal bone. 

Maxillary toothrow length.—Alveolar length of the upper 
cheek teeth. 

Nasal length/Occipitonasal length 
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Fic. 3. Three standards used for scoring skins for degree of develop- 
ment of lateral line. 


Mastoid breadth/Occipitonasal length 

Bullar length/Occipitonasal length 

Interparietal width/Mastoid breadth 

Maxillary toothrow length/Occipitonasal length 

Total length.—Recorded from skin tags. 

Tail length.—Recorded from skin tags. 

Hind foot length.—Recorded from skin tags. 

Tail length/Total length 

Hind foot length/Total length 

Background color.—Color of the dorsal underfur scored 
from 1 (light pinkish) to 5 (bright yellowish orange). 

Amount of dorsal black.—Color of the dorsal guard hairs 
scored from 1 (light) to 5 (dark). See figure 2. 

Postauricular patch.—Degree of development of post- 
auricular patches scored from 1 (none) to 5 (large). See 
figure 2. 
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Silkiness.—Softness of the pelage scored from 2 (silky) to 
4 (harsh). 

Lateral line.—Degree of development of lateral line scored 
from 1 (indistinct) to 3 (distinct). See figure 3. 


STATISTICAL TECHNIQUES 


After deleting specimens with missing characters, I sub- 
jected 57 P. merriami, 42 P. m. gilvus, and 37 P. flavus to a 
stepwise discriminant function analysis using a program 
(BMD07M) which originated at the Health Sciences Com- 
puting Facility, UCLA (Dixon, 1965). I treated the P. 
merriami and P. flavus as reference groups, and entered all 
42 P. m. gilvus as unknowns. A subsequent run using all three 
as distinct groups yielded substantially similar results. 

Explanations and examples of discriminant function anal- 
yses are available in Fisher (1936), Morrison (1967), 
Cooley and Lohnes (1962), and Sokal and Rohlf (1969). 
Basically, it is a technique for allocating unknown specimens 
to one of two or more previously recognized groups. It com- 
putes a new variable which is a linear function of the original 
n variables that maximizes the distance between the groups. 
Unknown specimens are scored for the new variable and 
placed in the most appropriate group. 

The analysis is done in stepwise fashion, adding or de- 
leting a single variable at each step. This allows one to 
determine which characters best separate the groups when 
used singly or in small combinations. In addition, within 
groups covariance and correlation matrices are generated, 
permitting one to determine the degree to which characters 
are correlated. 

The end product is a canonical analysis using all variables 
(except those which are deleted by the computer). A number 
of canonical variables equal to the number of characters used 
in the analysis is extracted. Character coefficients, which 
may be multiplied by their pooled within group standard 
deviation to show the contribution of each of the original 
variables to each of the new canonical variables, are given. 
In addition, each specimen is assigned a number corresponding 
to its generalized taxonomic distance ( Mahalinobis' distance) 


The Status of Perognathus merriami 185 


TABLE 3. Number of individuals classified as P. merriami or P. flavus 
after various steps in the discriminant function analysis 


Group merriami flavus 
Step number 1 merriami 53 4 
Variables used 1 flavus 6 31 
gilvus 35 7 
Step number 4 merriami 55 2 
Variables used 1—4 flavus 1 36 
gilvus 33 9 
Step number 8 merriami 56 1 
Variables used 1-8 flavus 0 ey 
gilvus Dn 15 
Step number 12 merriami 56 1 
Variables used 1-12 flavus 0 37 
gilvus 25 17 
Step number 16 merriami 56 1 
Variables used 1-16 flavus 0 37 
gilvus 24 18 
Step number 25 merriami 57 0 
Variables used 1-23 flavus 0 37 
gilvus gi 15 


from the mean of each reference group. The probability of 
a given specimen belonging to a given group is indicated for 
every specimen in all groups. 

Eigenvalues, cummulative proportion of total dispersion, and 
canonical correlations are calculated for each canonical vari- 
able. These values allow one to determine how much of the 
total available variability is accounted for by each successive 
canonical variable. 


RESULTS 


To demonstrate the difficulty encountered in attempting to 
separate these mice using morphological characters, Figure 4 
shows histograms for four of the better characters. These 
figures clearly show the overlap between P. merriami and 
P. flavus, and demonstrate the intermediate nature of P. m. 
gilvus. Table 2 is a list of the original variables in the order 
the computer entered them into the analysis. Table 3 shows 
the grouping results of the stepwise analysis at various steps. 
The characters contributing the most to the first canonical 
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TABLE 4. Standardized character coefficients and percent contributions 
to the first canonical variable 


Coefficient Percent 
Occipitonasal length 4.40 11.19 
Total length -3.83 9.74 
Interparietal width/Mastoid breadth 3.62 9.21 
Bullar length -3.60 9.16 
Interparietal width -3.56 9.05 
Bullar length/Occipitonasal length QU 8.16 
Hind foot length/Total length —3.07 7.81 
Nasal length -2.71 6.89 
Hind foot length 2.50 6.36 
Tail length/Total length -2.03 5.16 
Tail length 1.83 4.65 
Mastoid breadth 1.28 3.26 
Nasal length/Occipitonasal length 1.00 2.54 
Mastoid breadth/Occipitonasal length 0.45 1.14 
Lateral line 0.42 1.07 
Background color —0.39 0.99 
Postauricular patch 0.32 0.81 
Dorsal black —0.31 0.79 
Interparietal length 0.26 0.66 
Frontonasal length 0.24 0.61 
Maxillary toothrow length 0.22 0.56 
Silkiness —0.04 0.10 
Interorbita] constriction —0.03 0.08 


variable are not the same as those at the top of the list. Table 
4 gives the standardized coefficients and the percentage con- 
tribution of each of the original variables to the first canonical 
variable. 

Histograms of the discriminant scores of the first canonical 
variable for P. merriami, P. flavus, and P. m. gilvus are pre- 
sented in Figure 5. It can be seen that although complete 
separation of P. merriami and P. flavus can be achieved using 
multivariate techniques, once again P. m. gilvus bridges the 
gap. Figure 6 is a plot of the first and second canonical 


> 


Fic. 4. Histograms showing single character separation for four of 
the best characters. Perognathus flavus on top, P. m. gilvus in the 
middle, and P. merriami on the bottom. 
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P flavus 


n = 37 
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First Canonicol Variable 


Fic. 5. Histograms showing distribution of the first canonical variable. 
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Fic. 6. Plot of the first and second canonical variables. m = P. 
merriami, g = P. m. gilvus, f = P. flavus, * = group means, $ = overlap. 


variables, and includes all of the specimens used in the 
analyses. 


Discussion 


Table 2 and Figure 4 both illustrate the difficulties which 
faced taxonomists working with these animals in the past. 
Although differences can be seen between certain individuals, 
quantifying those differences and finding characters which are 
useful in allocating all individuals is extremely difficult. Table 
3 shows that the best single character, Bullar length/Occipito- 
nasal length, is not sufficient to allocate all of the individuals 
correctly, Using this character, four of the P. merriami and 
six of the P. flavus are incorrectly allocated, while 35 of the 
P. m. gilvus are allocated to P. merriami and seven to P. flavus. 
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Using the four characters shown in Table 2? in combination, 
only two P. merriami and one P. flavus are misallocated, and 
two additional P. m. gilvus are switched to P. flavus. All of 
the P. flavus and all but one of the P. merriami are correctly 
allocated using eight characters. At this step, 27 P. m. gilvus 
are assigned to P. merriami and 15 to P. flavus. This is the 
same as the final proportions for P. m. gilvus. At 16 and 20 
characters, the proportion of P. m. gilvus allocated to P. flavus 
continues upward through 17 to 18. The final step correctly 
allocates all of the reference samples, and places 27 P. m. 
gilvus in P. merriami and 15 in P. flavus. 

Of the reference P. merriami sample, one had a probability 
of 0.953 and one of 0.999 of being in the right group, and all 
of the rest were 1.000. For P. flavus, there were 2 with 0.953, 
1 with 0.993, 1 with 0.998, 2 with 0.999 and the rest were 
1.000. The 27 P. m. gilvus allocated to P. merriami had prob- 
abilities ranging from 0.953 to 1.000. Three P. m. gilvus had 
probabilities of 0.500 for both groups. 

The three intermediate specimens included one from 25 
mi. W of Tucumcari, New Mexico, one from Comstock, Texas, 
and one from near Kermit, Texas. Osgood (1900) considered 
the Comstock specimen to be an intermediate between typical 
P. merriami and P. m. gilvus. The type-specimen of P. m. 
gilvus was allocated to P. flavus with a probability of 0.731. 

The cumulative proportions of total dispersion for the first 
two canonical variables were 0.99861 and 0.99971, respectively, 
and corresponding eigenvalues were 8.21030 and 0.00906. This 
indicates that the data illustrated in Figures 5 and 6 account 
for virtually all of the variability within these samples. 

Figure 6 probably illustrates quite well the actual relation- 
ships between these three groups. Perognathus m. gilvus 
is a geographically intermediate population between P. flavus 
to the west and P. merriami to the east. I suggest that the 
overlap shown in Figure 6 is probably indicative of the popu- 
lation now known as P. m. gilvus interbreeding with P. flavus 
in Eastern New Mexico and Western Texas, and with P. 
merriami in south-central Texas. Data are at present in- 
sufficient to determine whether the zone of overlap between 
P. flavus and P. m. gilvus is a narrow one limited to Eastern 
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New Mexico and adjacent portions of Western Texas or 
whether the integradation is more broadly clinal over the 
whole range of P. m. gilvus. 

Throughout this analysis, P. m. gilvus has been demonstrated 
to be closer to P. merriami than to P. flavus. This is not 
surprising, since the majority of the specimens are from areas 
that are geographically closer to P. merriami than to P. flavus. 
This probably also explains the early allocations of this popu- 
lation to P. merriami, rather than to P. flavus, even though the 
close relationships of all three were noted (Osgood, 1900). 

In conclusion then, it seems best to regard P. flavus and 
P. merriami as conspecific. The three populations are cer- 
tainly distinct enough to retain subspecific status with the 
easternmost population known as Perognathus flavus merriami, 
and the intermediate population as Perognathus flavus gilvus. 
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